Study Area: GoMex

library(sp)
library(rgdal)
library(raster)
shore = shapefile('downloads/GSHHS_f_GOM/GSHHS_f_GOM.shp')
#plot(shore)

Raw Impacts

download

library(dataone) # install.packages("dataone")

options(dataone_token = readLines('.dataone_token'))
am = AuthenticationManager()
getTokenInfo(am)
cn = CNode("PROD")
q = query(cn, solrQuery=list(q="id:doi*", fq="id:doi\\:10.5063/F1S180FS"))
m = q[[1]]

# download tifs
for (doc in m$documents){ # doc = m$documents[[2]]
  url = sprintf('https://knb.ecoinformatics.org/knb/d1/mn/v2/object/%s', doc)
  zip = sprintf('downloads/%s.zip', doc)
  download.file(url, zip)  
  unzip(zip, exdir='downloads')
}

crop

# project shoreline from geographic to Mollweide like raster
r_1 = raster(list.files('downloads','.*\\.tif$', full.names = T)[1])
shore_mol = spTransform(shore, crs(r_1))

for (tif in list.files('downloads','.*\\.tif$', full.names = T)){ # tif = list.files('downloads','.*\\.tif$', full.names = T)[1]
  r_g = raster(tif)
  r_mol = crop(r_g, extent(shore_mol))
  writeRaster(r_mol, sprintf('tif_mol/%s', basename(tif)))
}

plot

# project shoreline from geographic to Mollweide like raster
r_1 = raster(list.files('tif_mol','.*\\.tif$', full.names = T)[1])
shore_mol = spTransform(shore, crs(r_1))

# loop through tifs
for (tif in list.files('tif_mol','.*\\.tif$', full.names = T)){ # tif = list.files('tif_mol','.*\\.tif$', full.names = T)[1]
  r = raster(tif)
  plot(
    r, main = tools::file_path_sans_ext(basename(tif)),
    col=rev(colorRampPalette(RColorBrewer::brewer.pal(n=9,'Spectral'))(256)))
  plot(shore_mol, add=T)
}